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Abstract 

(N 

^ I SFITTER is a new analysis tool to determine supersymmetric model parameters 

^ ■ from collider measurements. Using the set of supersymmetric mass measurements 

^ , at the LC and at the LHC we show how both colliders probe different sectors of the 

tJ- I MSSM Lagrangian. This observation is a strong motivation to move from a parameter 

fit assuming a certain model to the unconstrained weak-scale MSSM Lagrangian. We 



o 

^ . argue how the technical challenges can be dealt with in a combined fit/ grid approach 

with full correlations. 

&: 1 Introduction 

While the Standard Model describes all available high energy physics experiments, it 
^ . still has to be regarded as an effective theory, valid at the weak scale. New physics are 
^ ! expected to appear at the TeV energy scale. The supersymmetric extension ( [1]) of the 
Standard Model is a well motivated extension providing us with a description of physics 
that can be extended consistently up to the unification scale. 

If supersymmetry or any other high-scale extension of the Standard Model is discov- 
ered, it will be crucial to determine its fundamental high-scale parameters from weak- 
scale measurements [2]. The LHC and future Linear Colliders will provide us with a 
wealth of measurements [3, 4], which due to their complexity require proper treatment 
to unravel the corresponding high-scale physics. Even in the general weak-scale mini- 
mal supersymmetric extension of the standard model (MSSM [5]) without any unifica- 
tion or SUSY breaking assumptions some of the measurements of masses and couplings 
are not independent measurements; moreover, linking supersymmetric particle masses 
to weak-scale SUSY parameters involves non-trivial mixing to mass eigenstates in essen- 
tially every sector of the theory. On top of that, for example in gravity mediated SUSY 
breaking scenarios (mSUGRA/cMSSM) a given weak-scale SUSY parameter will always 
be sensitive to several high-scale parameters which contribute through renormalization 
group running. Therefore, a fit of the model parameters using all experimental infor- 
mation available will lead to the best sensitivity and make the most efficient use of the 
information available. 
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Table 1: Errors for the mass determination in SPSla, taken from [17]. Shown are the 
nominal parameter values and the error for the LHC alone, the LC alone, and a combined 
LHC+LC analysis. All values are given in GeV. 

In a fit, the allowed parameter space might not be sampled completely. To avoid 
boundaries imposed by non-physical parameter points, which can confine the fit to a 
'wrong' parameter region, combining the fit with an initial evaluation of a multi-dimens- 
ional grid is the optimal approach. 

In the general MSSM the weak-scale parameters can vastly outnumber the collider 
measurements, so that a complete parameter fit is not possible and one has to limit one- 
self to a consistent subset of parameters. In SFITTER both grid and fit are realised and 
can be combined. This way, one can ultimately elim.inate all dependence on the starting 
point of the parameter determination. SFITTER also includes a general correlation matrix 
and the option to exclude parameters of the model from the fit/ grid by fixing them to 
a value. Additionally, SFITTER includes the option to apply a Gaussian smearing to all 
observables before they enter the fit/grid in order to simulate realistically experimiental 
measurements. In this preliminary study, however, correlations and systematic uncer- 
tainties are neglected and the central values are used for the measurements. 

Currently, SFITTER uses the predictions for the supersymmetric masses provided by 
SUSPECT [6], but the conventions of the SUSY Les Houches accord [7] allow us to in- 
terface other programs. The branching ratios and e'^e" production cross sections are 
provided by MSMlib [8], which has been used extensively at LEP and cross checked 
with Ref. [9]. The next-to-leading order hadron collider cross sections are computed 
using PROSPINO [10]. The fitting program uses the MINUIT package [11]. The deter- 
mination of includes a general correlation matrix between measurements. In its next 
version SFITTER will be interfaced with the improved branching fraction determination 
of SDECAY [12], as well as alternative renormalization group codes like SoftSUSY [14], 
ISAJET [15] or SPHENO [16]. 
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Table 2: Summary of the mSUGRA fits in SPSla: true values, starting values, fit values 
and absolute errors from the fit. As in SPSla we fix /x > 0. The mass values of the fits are 
based on can be found in Tab.lll 

2 mSUGRA/cMSSM Parameter Determination 

Assuming that SUSY breaking is mediated by gravitational interactions (mSUGRA/ cMSSM) 
we fit four universal high-scale parameters to a toy set of collider measurements: the uni- 
versal scalar and gaugino masses, mo, mi/2, the trilinear coupling Ao and the ratio of the 
Higgs vacuum expectation values, tan /3. The sign of the Higgsino mass parameter /i is a 
discrete parameter and therefore fixed. In contrast to an earlier study [13] we assume the 
set of mass measurement at the LHC and at the LC, shown in Tab. [H The central value 
for our assumed data set corresponds to the SUSY parameter point SPSla [18], as com- 
puted by SUSPECT. As mentioned in the introduction correlations, systematic errors and 
theoretical errors are neglected. As the central (true) values are used as measurements in 
order to study the errors on the determination of the parameters, the values are not 
meaningful and therefore are not quoted. 

The starting points for the mSUGRA parameters are fixed to the mean of the lower 
and upper limit (typically 1 TeV/c^) of the allowed parameter range, i.e. they are not 
necessarily close to the true SPSla values. The result of the fit is shown in Tab. |2| All 
true parameter values are reconstructed well within the quoted errors, in spite of starting 
values relatively far away. The measurements of mo and mi/2 are very precise, while the 
sensitivity of the masses on tan P and Aq is significantly weaker. The results for the LHC 
alone are generally an order of magnitude less precise than those for the LC, and this 
qualitative difference is expected to become even more pronounced once we properly 
include systematical errors. 

Because the data set is fit assuming mSUGRA as a unification scenario the absence of 
measurements of most of the strongly interacting particles, in particular the gluino, does 
not have a strong impact on the precision of the LC determination. Therefore the results 
for the combined measurements LHC+LC show only a small improvement. 

Assuming an uncorrelated data set, the correlations between the different high-scale 
SUSY parameters which we obtain from the fit are given in Tab. |3| We can understand 
the correlation matrix step by step [19]: first, the universal gaugino mass mi/2 can be 
extracted very precisely from the physical gaugino masses. The determination of the 
universal scalar mass mo is dominated by the weak-scale scalar particle spectrum, but 
in particular the squark masses are also strongly dependent on the universal gaugino 
mass, because of mixing effects in the renormalization group running. Hence, a strong 
correlation between the mo and mi/2 occurs. The universal trilinear coupling Aq can be 
measured through the third generation weak-scale mass parameters A^^t.T- However, the 

t T- which appear for example in the off-diagonal elements of the scalar mass matrices, 
also depend on mo and mi/2/ so that Aq is strongly correlated with mo and mi/2- At this 
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Table 3: The (symmetric) correlation matrix for the mSUGRA fit given in TabJUwith data 
set LHC+LC. 

point one should stress that the determination of Aq is likely to be dominated by At as 
it appears in the calculation of the lightest Higgs mass ruh- After taking into account 
the current theoretical error of 3 GeV on rrih [20] we expect the determination of Aq to 
suffer significantly The experimental errors therefore can be considered a call for an 
improvement of the theoretical error. 

In general, tan P can be determined in three sectors of the supersymmetric spectrum: 
all four Higgs masses, and for large values of niA in particular the light CP even Higgs 
mass ruh depend on tan/?. The mixing between gauginos and Higgsinos in the neu- 
tralino/chargino sector is governed by tan/5. Finally, the stop mixing is governed by 
fi/ tan P, while the sbottom and stau mixing depends on /i tan /?. The correlation of tan P 
with the other model parameters reflects the relative impact of these three sectors. In an 
earlier analysis we assumed a uniform error of 0.5% on all mass measurements [13] and 
saw that in this case tan P is determined through stau mixing, which in turn means that it 
shows very little correlation with mi/2- 

For the more realistic scenario in Tab. [l] the outcome is the following: the relative 
errors for the light Higgs mass and for the light neutralino masses at the LC are tiny. The 
relevant parameter in the Higgs sector is the light stop mass, which is governed by mi/2; 
similarly the gaugino mass mi/2 which fixes the light neutralino and chargino masses 
does not depend strongly on tan p. The slepton sector introduces a strong correlation 
between niQ and mi/2. The resulting correlation matrix is shown in in Tab.|3l The results 
obtained with SFITTER are in agreement with expectation. 

3 General MSSM Parameter Determination 

In this study, the unconstrained weak-scale MSSM is described by 24 parameters in ad- 
dition to the standard model parameters. The parameters are listed in Tab. HI tan P as 
in mSUGRA, plus three soft SUSY breaking gaugino masses Mj, the Higgsino mass pa- 
rameter /i, the pseudoscalar Higgs mass m^, the soft SUSY breaking masses for the right 
sfermions, Mj^, the corresponding masses for the left doublet sfermions, Mj^ and finally 
the trilinear couplings of the third generation sfermions At^b,T- 

3.1 Toy model with all masses 

For testing purposes, we first consider a toy data set which includes all supersymmetric 
particle masses. The universal error on all mass measurements is set to 0.5%. 

In any MSSM spectrum, in first approximation, the parameters Mi, M2, /x and tan P 
determine the neutralino and chargino masses and couplings. We exploit this feature to 
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Table 4: Result for the general MSSM parameter determination in SPSla using the toy 
sample of all MSSM particle masses with a universal error of 0.5%. Shown are the nominal 
parameter values, the result after the grid and the final result. All masses are given in GeV. 

illustrate the option to use a grid before starting the fit. The starting values of the param- 
eters other than Mi, M2, and tan/3 are set to their nominal values, this study is thus 
less general than the one of mSUGRA. The is then minimized on a grid using the six 
chargino and neutralino masses as measurements to determine the four parameters Mi, 
M2, fi and tan p. The step size of the grid is 10 for tan P and 100 GeV for the mass pa- 
rameters. After the minimization, the four parameters obtained from grid minimization 
are fixed and all remaining parameters are fitted. In a final run all model parameters are 
released and fitted. The results after the grid (including the complementary fit), after the 
final fit and the nominal values are shown in Tab. HI The smearing option has not been 
applied. However, the errors on the fitted values (once the fit converges) should not be 
sensitive to these shortcomings. 

The final fit values indeed converges to the correct central values within its error. The 
central values of the fit are in good agreement with generated values, except for the trilin- 
ear coupling Ab^r- The problem is using only mass measurements to determine the three 
entries in a (symmetric) scalar mass matrix: in the light slepton sector there are three 
masses, left and right scalars plus the sneutrino, so the system is in principle calculable. 
In the third generation squark sector we have three independent diagonal entries per 
generation and two off-diagonal entries. But the number of mass measurements is only 
four, therefore the system is underdetermined in first order. The off-diagonal entry in the 
mass matrix for down type scalars includes a term A^^r and an additional term n tan p. 
Even for very moderate values of tan P the extraction of Af,^ requires precise knowledge 
of tan p. The use of branching ratios and cross section measurements (with polarised 
beams) which carry information about the scalar mixing angles should significantly im- 
prove the determination of At^b,T- 

3.2 Toy model with LHC+LC mass measurements 

In the study of the three data sets LHC, LC, and LHC+LC in the MSSM, a fit was per- 
formed for the data sets LHC and LC, whereas for LHC+LC additionally the GRID was 
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Table 5: Result for the general MSSM parameter determination in SPSla using the mass 
measurements given in Tab.[ll Shown are the nominal parameter values and the result 
after fits to the different data sets. All masses are given in GeV. 

used for Mi, M2, /i and tan (3 with the five chargino and neutralino masses. The starting 
points were chosen to be the true values (with the exception of the parameters used in 
the grid). In order to obtain a solvable system, for the LHC data set m^, Mi^, Mf^, 
were fixed. For the LC data set the first and second generation squark soft SUSY breaking 
masses, the gluino mass M3, Mj^ and Ai, were fixed. These parameters were chosen on 
the basis of the measurements available in Tab. [H The values to which these parameters 
were fixed is not expected to influence the final result of the fit. The results for the two 
data sets are shown in Tab.|51 

Note that the general rule that the LHC is not sensitive to weakly interacting particle 
masses is not entirely true: while the LHC has the advantage of measuring the squark 
and gluino masses, the first and second generation slepton mass parameters are also de- 
termined with a precision of the order of percent. The results in Tab. |5l show that the LHC 
alone is well capable of determining for example all gaugino mass parameters as well as 
most of the scalar mass parameters. 

The situation at the LC is slightly different. Only marginal information on the squark 
sector available at the LC. The measurement of At from the Higgs sector should be taken 
with a grain of salt (theoretical error on the lightest Higgs mass). Adding the stau mixing 
angle to the set of LC measurements will improve the determination of A^.. However, the 
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measurements of the parameters, in particular slepton and gaugino parameters are far 
more precise than at the LHC. 

For the LHC+LC data set, a sufficient number of mass measurements is available, so 
that no parameters need to be fixed. The superiority of the combination of the measure- 
ments at the two colliders is obvious from this observation and from Tab.|5| The LHC con- 
tributes to reduce the error in the weak sector (A/2) and the LC in the strongly interacting 
sector (third generation squarks). Even more important: of 13 parameters undetermined 
by either the LHC or the LC, 11 are determined with good precision in the combination. 
For At, Ab, we expect an improvement with the use of branching ratios and cross section 
measurements. 

A complete measurement of all parameters at the weak scale is particularly important 
if one wants to probe unification scenarios which link subsectors of the parameter space 
which are independent at the weak scale. An advanced tool like SFITTER can extract the 
information to probe supersymmetry breaking scenarios from any set of measurements, 
provided the set is sufficient to overconstrain the model parameters. 

4 Conclusions 

SFITTER is a new program to determine supersymmetric parameters from measurements. 
The parameters can be extracted either using a fit, a multi-dimensional grid, or a combi- 
nation of the two. Correlations between measurements can be specified and are taken 
into account. While it is relatively easy to fit a fixed model with very few parameters 
for example at a high scale to a set of collider measurements, the determination of the 
complete set of weak-scale MSSM model parameters requires this more advanced tool. A 
mSUGRA inspired fit does not include the full complexity and power of the combined 
LHC and LC data compared to the measurements at either collider alone. The results 
from SFITTER in the MSSM with the three data sets show that only the combination of 
measurements of both the LHC and the LC offers a complete picture of the MSSM model 
parameters in a reasonably model independent framework. 
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